The preliminary lattice QCD calculation of k meson decay width 
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We present a direct lattice QCD calculation of the k meson decay width with the s-wave scattering 
phase shift for the isospin 1 = 1/2 pion-kaon (tt/sT) system. We employ a special finite size formula, 
which is the extension of the Rummukainen-Gottlieb formula for the nK system in the moving frame, 
to calculate the scattering phase, which indicates a resonance around k meson mass. Through the 
effective range formula, we extract the effective k. — iiK coupling constant Qk-wK ~ 4.54(76) GeV 
and decay width F — 293 ± 101 MeV. Our simulations are done with the MILC gauge configurations 
with Nf =2+1 flavors of the "Asqtad" improved staggered dynamical sea quarks on a 16"^ x 48 
lattice at (m-^ + mK)/mK ~ 0.8 and lattice spacing a ~ 0.15 fm. 
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I. INTRODUCTION 

It is well-known that kappa meson is a resonance, 
which is a state with a considerable width under strong 
interactions, and decay via strong interaction. In 2010, 
Particle Data Group (PDG) [l| lists the meson i^i; (800), 
which is usually called k meson {J^ = O"*", / = 1/2), 
with a mass (676 ± 40 MeV) and a broad width (548 ± 24 
MeV). Some recent experimental analyses |2l-[l2| and 
phenomenological analysis of nK scattering 
strongly show its existence. The Dalitz plot analysis 
of Fermilab experiment E791 Q yield its decay width 
about 410 ± 43 ± 87 MeV. Moreover, the k meson has 
been extensively studied with BES data 0, Q , where the 
evidence for its existence is very clear, and the most re- 
cent analyses based on events collected by BESII gives 
its mass about 849 ± 77 MeV, and decay width about 
512 ± 80 MeV. 

Although the direct determination of k resonance pa- 
rameters from QCD is extremely hard since the calcula- 
tion of resonance masses and decay widths is essentially 
a nonperturbative problem. Several research groups have 
undertaken theoretical efforts to study the k meson and 
its resonance parameters However, all the exper- 

imental and theoretical analyses give a little bit different 
resonance mass and decay width for k meson. Therefore, 
there is not close to the consensus yet on its resonance 
parameters. 

The feasible way to study k resonance nonperturba- 
tively from first principles is with the help of lattice 
QCD. Until now, no direct lattice QCD study about k 
resonance have been reported, possibly because the reli- 
able calculation of the rectangular diagram are difficult, 
and there arc not suitable theoretical formula available 
to represent the ttK system in the moving frame. In- 
spired by J. Ncbreda and J. Pelaez's theoretical studies 
on K resonance [l^ and our previous works about precise 
extraction of k mass [131, and reliable extraction of ttK 
scattering length |2l| at / = 1/2 channel, in this work, 
we further explore its resonance parameters directly from 
lattice QCD. 



In this paper, we estimate the k meson decay width 
by calculating the s-wave scattering phase shift for the 
7 = 1/2 ttK system. We discuss them both in center of 
mass (CM) frame and the moving frame. The simulation 
are done with MILC lattice ensemble [12, [1^ • The meson 
masses determined in [13| gave (m^-l-mif)/™^ = 0.8 [13], 
and the lattice parameters determined by MILC collab- 
oration are: the lattice extent L « 2.5 fm and the lattice 
space inverse 1/a = 1.358 GeV [H, |11|. The finite size 
formula (23 - [2^ is employed to nK system in the center- 
of-mass frame, and we utilize a special finite size formula, 
which is the extension of the Rummukainen-Gottlieb for- 
mula for the ttK system in the moving frame [s^ . to 
compute the scattering phase shift. 



II. METHODS 

A. The effective range formula 

In the ttK system, the relativistic Breit-Wigner for- 
mula (RBWF) for the elastic s-wave scattering phase i5o 
in the resonance region with a center-of-mass energy to^ 
and a decay width F^ can be conveniently expressed as 



tanJn = 



/5Ffl(s) 



(1) 



where s = E^^^j is the invariant mass of the ttK sys- 
tem, and EcM is its center-of-mass energy. The quantum 
numbers of k resonance is I{J^) = ^(0^) and decays into 
one pion and one kaon in the s-wave, which can be han- 
dle on the lattice. The decay width Fh(s) can be written 
in terms of the coupling constant f^KTrif [lal, 



ri?(s) 



8tt s ' 



(2) 



where 



P = 7ri= v^[s - {rriyr - mif)2] [s - (m^ -f- 771^)2] , (3) 
2Vs 

is the center-of-mass momentum of the pion or kaon. 
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With a combination of Eqs. ([T|) and a description 
of the scattering phase in the s-wave as a function of the 
invariant mass ^/s is provided by so-cahed the effective 
range formula (ERF) in the elastic region, 



tan. So 



Stt \fs(m}j^ — s) 



(4) 



This equation allows us a fit or seeking for the desired 
parameter gK^K, which is the effective k — > ttK coupling 
constant, and the resonance position m^. 

The K decay width r„ can then be calculated through. 



representation r = T+. The finite -size formula connect- 
ing the energy E or E to the s-wavc t:K scattering phase 
shift 60 is [3[]2i] 



tan (5o(fc) 



^00 (1; 9^)' 

where the zeta function is formally denoted by 

1 



g2) 



(7) 



(8) 



r.=r^(s) 
1 



9kttK Pii 

Stt ml ' 



(5) 



Pk 



2m. 



y/[inl - {rrij, - mKY][m1 - (m^ -|- m^)^ 



Therefore, equations ([4]) and ([5]) allow us to archive 
and by way of the dependence of 5q on s. 



B. Finite- volume method 

1. Center of mass frame 

In the center-of-mass frame, if we doesn't consider the 
interaction between pion and kaon, the possible energy 
eigenvalues of ttK system reads 



E = ^ml + p2 + y' m|. + , 

where p = |p|, p = (27r/L)n, and n g Z'^. In a typi- 
cal lattice study, the energy for the n ^ is larger than 
kappa mass TOk- For instance, on our chosen gauge con- 
figurations, the lowest energy for the n 7^ evaluated 
from TOtt, niK and to„ [2^, is i? = 1.35 x m^, which is 
not eligible to study the k meson decay. So, we must 
consider the n = case. On our chosen gauge config- 
urations the invariant mass of free pion and kaon takes 
y/s = 0.8 X itIk, which is away from kappa mass 77i„, but 
closer to itik than those with the n 7^ 0. 

With the consideration of the interaction of ttK sys- 
tem. The energy eigenvalues are moved from E to E, 
and the energy eigenvalue for ttK system is. 



E 



ml + fc2 + y^m|. + 



27r 
1 



where g e M. From this equation, we obtain the scatter- 
ing momentum k 



1 ^2 

k=-=y[E ~ {m^-mK)'^][E -{m^+mxy] . (6) 
2E 

In the center-of-mass frame these energy eigenstates 
transform under the cubic group 0(3) in the irreducible 



2. Moving frame 

In order to make the energy of ttK system is more 
close to the k mass m„, we consider the moving frame 
(MF) [l^l- Using a moving frame with total momentum 
P = (27r/L)d, d € Z'^, the energy eigenvalues for non- 
interacting irK system are 



Emf = Jml+pf + Jmj^+p^, 



(9) 



where pi = |pi|, P2 = IP2I, and pi, P2 denote the three- 
momenta of the pion and kaon meson, respectively, which 
satisfy the periodic boundary condition. 

Pi = -^n, , Hi e Z'^ , 

and the relation pi + p2 = P . 

Using the standard Lorcntz transformation with a 
boost factor 7 = l/Vl — v^, here v = P/Emf- The 
EcM can be obtained through 



EcAi--"f-^EMF = \^rn^+P*^+^Jml+p*-2, (10) 

where, in the center-of-mass frame, the total center-of- 
mass momentum disappears. 



P 



Pi = P2 



here and later we denote the center-of-mass momenta 
with an asterisk (*). The boost factor acts in the direc- 
tion of the velocity v, here we use the shorthand notation 

7P = 7P||+p±, 7"^p = 7"^P|] + P± , (11) 

where p|| and pj^ are the components of p parallel and 
perpendicular to the center-of-mass velocity, respectively. 



P|l = ^^' P±=P-P|| 



The p* are quantized to the values [35 



. 27r 
P =^r, 



rePd 



(12) 



(13) 
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where the set Pd is 



Correlation matrix 



Pd 



2'^ 



n e Z'' 



(14) 

Here we only consider the dominant low energy states 
in the moving frame: the pion with zero momentum, and 
kaon with the momentum p = (2?:/ L)e^ (d = 63) and the 
K meson with the momentum P = p. For our case, the 
invariant mass of ttK system has ^/s = 0.93 x m^, which 
is much closer to k mass than these in the center- 
of-mass frame. So, we only investigate this case. The 
invariant mass ^/s was calculated by ^/s = \/Ej^jp- 

In the interacting case, the Ecm is provided by 



P2. 



27r 



EcAi = V^l + k^ + J m|. + , k=—q. (15) 



L 



where q G M. From this equation, we gain 



k=^^J[E^-{m^-mK)^][E''-im^ + mKf] . (16) 
2E 



An immediate consequence of Eq. (|16p is that the physi- 
cal kinematics for the K-mcson decay satisfy 



Ecm > m^r + rriK, 



(17) 



namely, k can decay only if its mass exceeds the sum of 
the masses of its decay products (namely, tt and K). We 
can rewrite equqtion (jl6p to an elegant form for later use 



K 



(18) 



The energy eigenstates transform under the tetragonal 
group Civ the irreducible representations of Ai and E 
are relevant for the s-wave ttK scattering states in a cubic 
box. In this work, we only consider the energies linked 
with the Ai sector. The hadron interaction moves the 
energy from E to E, and the energies E are associated 
to the ttK scattering phase shift Sq in a torus through the 
ttK system's Rummukainen- Gottlieb formula [30j . which 
is an extension of the Liischer formula [2J] to the moving 
frame 1321. 



tan(5o(fc) 



2o1,(i;<z^) 



(19) 



where the modified zeta function is formally defined by 

1 



(20) 



rePd 



and the set Pd is denoted in Eq. p4| . Equation pOj) is 
first provided in Ref. [Hj for generic two-particles system, 
and we further confirmed and rigorously proved it [3^ . 
The k is the momentum defined from the invariant mass 
^/s as ^/s = \fW- 



■rat. 



The calculation 



method of Zqq(1;(7^) is full discussed in Appendix El 



In order to calculate the two energy eigenstates, 
namely, E^ (n ~ 1,2), we construct a matrix of the time 
correlation function, 



C{t) 



(0|Ol^.(0O.x(0)|0) (0|Ol,^(t)O«(0)|0) 
(0|Ot(t)O.A-(0)|0) (0|Ot(t)O.(0)|0) 



(21) 

where (0 is an interpolating operator for the k meson 
with the definite momentum, and OTrxit) is an inter- 
polating operator for the ttK system with specified mo- 
mentum. The interpolating operators 0^(1) and 0.„K{t) 
employed here is exactly the same as in our previous stud- 
ies [20, [111- To make this paper self-contained, we will 
provide all the necessary definitions below. 



1. ttK sector 

Here we briefly review the required formulae to com- 
pute the scattering phase for the ttK system in a torus 
using the original derivation and notation in Refs. (ssl - 
[37| . Let us study the nK scattering system of one 
Nambu-Goldstone pion and one Nambu-Goldstone kaon 
in the Asqtad-improved staggered dynamical fermion for- 
malism. Using the operators ©^(xi), C'7r(a;3) for pions 
at points a;i,X3, and the operators Ok{x2),Ok{x4) for 
kaons at points a;2,a;4, respectively, with the pion and 
kaon interpolating fleld operators defined by 

0,+ (x,i) = -d(x,t)75u(x,t) 

0^0 (x,i) = -^[?l(x,i)75w(x,t) - d(x,t)75d(x,i)]. 



Ox«(x,t) = s(x,t)75(i(x,t), 
Ok+ (x, t) = s(x, t)75M(x, t). 



(22) 



We then represent ttK four-point functions as 

C^k(x4,X3,X2,Xi) = {O K{xi)0T:{x-i)0^Ki^2)0\{xi)) . 

In this paper we only consider one moving frame, which is 
the pion with zero momentum, and kaon with momentum 
p = (27r/L)e3 (namely, d = 63) and total momentum 
P = p. The operator which creates a single kaon with 
non-zero three momentum k from the vacuum is obtained 
by Fourier transform: ©^(k, t) = ©^(x, t)e^^'^ . 

After summing over the spatial coordinates xi , X2 , X3 
and X4, we get the ttK four-point function 

CttK (P, ^4, ^3, ^2, il) = 

e*P-(x4-x.)c.^^(^^^ X3, X2, xi) , (23) 

Xi ,X2 ,X3 ,X4 

in the momentum p state, where xi = (xi,ti), X2 = 
(x2,t2); x^ = (x3,t3), and x^ = (x4,t4), and t stands for 
the time difference, namely, t = — ti. To prevent the 
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(a) D (b) C fc) R 



FIG. 1: Diagrams contributing to ttA' four-point functions. 
Short bars stand for the wall sources. Open circles are sinks 
for local pion or kaon operators. The thicker lines represent 
the strange quark lines. 

complicated Fierz rearrangement of the quark lines [36j . 
we select ti = 0,t2 = 1,^3 = t, and <4 = t + 1. We 
construct ttK operators for the 7=1/2 channel as 

0ilHp,t) = i=|V27r+(t)i^"(p,i + l)- 

7T'>{t)K+{p,t + l)}, (24) 

where p is total momentum of ttK system. This 7:K 
operator has (7,/^) = (1/2,1/2). 



Supposing that the masses of the u and d quarks are 
degenerate, only three diagrams contribute to the ttK 
scattering amplitudes [131 . These quark line diagrams are 
shown in Fig.[TJ which are labeled as direct diagram (D), 
crossed diagram (C) and rectangular diagram (R), re- 
spectively. The direct and crossed diagrams can be read- 
ily evaluated [s^. However, the rectangular diagram (R) 
needs another quark propagator linking the time slices t^, 
and ^4, which make it hard to compute (36| . 

Sasaki ct al. tackle this puzzle through the technique 
with a fixed kaon sink operato r 1381 . We successfully 
extended the technique in Refs. [SJ, [Sg] to the ttK scat- 
tering at the 7=1/2 channel with zero momentum [2l| . 
here we further consider the case with the non-zero mo- 
mentum. To be specific, each propagator corresponding 
to a moving wall source is given by [sj, [3^ 

Y,Dn',n"Gtin")=J2^"',i^,t), 0<t<T-l. (25) 

n" X 

For the non-zero momentum, we used a up quark 
source with 1, and a strange quark source with e^^ '' on 
each site for the pion and kaon creation operator, respec- 
tively. 7?, C, and R, arc shown in Fig. [1] we can express 
them in terms of the quark propagators G, namely. 



Cf^(p,i4,t3,t2,tl) = e'P-"^(RcTr[Gl^(x3,t3)Gt,(x3,t3)G'Ux4,t4)Gt,(x4,i4)] 

X3 X4 

C^j,{p,h,h,t2,h) = e*P-"'(ReTr[Gj^(x3,t3)Gt,(x3,t3)G'Ux4,i4)Gt,(x4,t4)] 

X3 X4 

C^Kip,U,ts,t2,h) = e*P'"'(RcTr[Gj^(x2,t2)Gt,(x2,t2)Gl^(x3,t3)Gt,(x3,t3)] 

I 



(26) 



All three diagrams in Figure [T] are needed to study ttK 
scattering in the 7=1/2 channel. If assuming that u 
and d quarks have the same mass, the ttK correlation 



2. K sector 



In our previous work (20[, we make a detailed proce- 



function for 7 = 1/2 channel can be expressed as the dure to measure (0|Kt(t)K(0)|0). For the light u quark 



combinations of these three diagrams [3 



Dirac operator M„ and the s quark Dirac operator Mg, 
we obtain k correlator [2G| 

1 3 

C^Kit) EE {O^Kit)\O^Km =D + -C--R, (27) ^ ^(-l)-e'P-"(Tr[M-i(x,<;0,0)M,-i'(x,i;0,0)]), 

(29) 

with given momentum p. 

For staggered quarks, the meson propagators have the 
generic single-particle form, 

Cit) = ^ A.e-™'* + ^ Aa-l)*e-™-* + {t-^Nt~ t), 

i i 

(30) 

where 0,^(0, t) and Gif(l,i -I- 1) are the pion and kaon where the oscillating terms correspond to a particle with 



where the operator O^x denoted in Eq. creates a 
ttK state with total isospin 1/2. 

In our concrete calculation we also calculate the ratios 



two-point functions, respectively. 



opposite parity. For the k meson correlator, we consider 
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only one mass with each parity in Eq. (|30l) . and the k 
correlator was fit to 

= 6„e-™''* + &K,(-)*e-*^--* + (t ^ 7Vi-t), (31) 

where the ^-^d two overlap factors. 

3. Off-diagonal sector 

To avoid the complicated Ficrz rearrangement of the 
quark lines [s^, we select ti = 0,t2 = 1, and ts = < for 
ttK — >■ K three-point function, and choose ti = 0,^2 = t, 
and tz = t + I for k — >■ tt^sT three-point function. 




FIG. 2: Diagrams contributing to nK — >■ and k — > Trie" 
three-point functions, (a) Quark contractions of ttK — > k 
and (b) Quark contractions of k — vrif. Short bars show the 
wall sources. Open circles are sinks for local pion operator. 
The thicker lines indicate the strange quark lines. 

The quark line diagrams contributing to /t — >■ ttK and 
ttK —7- K three-point function arc plotted in Fig. [H^a) 
and Figl2fb), respectively. For the nonzero momentum, 
we used a u quark source with 1, and a strange quark 
source with e^^ '' on each site for kaon creation operator. 
We can express t:K — > k three-point function in terms of 
the quark propagators G as 

J2 e^^-'^^ (Tr[75Gt, i^3,t3)Gl (xg, h)^5Gl (xi , ti)] ,(32) 

where trace is over the color index; the 75 matrix are 
utilized as an interpolating field for the pseudoscalars. 



D. Extraction of energies 

To reliably obtain two lowest energy eigenstates, we 
employ the variational method [13] and construct a ratio 
of the correlation function matrices as 

M{t,tR)^C{t)C-\tn), (33) 

with some reference time slice tu (40| . Two lowest energy 
levels can be obtained by a correct fit to two eigenvalues 
^n{t,tii) (n = 1,2) of the matrix M{t,tji). Since we 



work on the staggered fermion, \n{t,tji) [n = 1,2) are 
behaved as [2lj 

Ht,tR) = A„cosh{~E„{t-T/2)) 

+ (-l)*i?„cosh(-<(t-r/2)), (34) 

which explicitly contains an oscillating term, for a large 
t, which mean <C ti?, < i ^ T/2 to suppress the excited 
states and the wrap-around contaminations [13]. Here 
we suppose the non-dcgencratc eigenvalues Ai(t,ii?) > 

X2{t,tR). 

III. LATTICE CALCULATION 

We used the MILC lattice ensemble with the Nf = 2+1 
dynamical flavors of the Asqtad-improved staggered dy- 
namical fermions, see the more details in Refs. [2^ l23j. 
We measured ttK four-point functions on the 0.15 fm 
MILC lattice ensemble of 360 16'^ x 48 gauge config- 
urations with bare quark masses arriud = 0-0097 and 
aiTis = 0.0484 and bare gauge coupling 10/^^ = 6.572, 
which has a physical volume approximately 2.5 ftn, and 
the inverse lattice spacing 0"^ = 1.358l^^ GeV 0,[2ll. 
The mass of the dynamical strange quark is near to its 
physical value, and the masses of the u and d quarks have 
the same mass. 

We use the standard conjugate gradient method to get 
the required matrix element of the inverse fermion ma- 
trix. The correlator Cii(t) is calculated by 

Ciiit) = ^5]((^i^)(t + t.)(^i^)Ut.)). (35) 

After averaging the correlator over all 48 possible values, 
the statistics are greatly improved. 

We calculate off-diagonal correlator C2i{t) through 

C2i{t) = {nit)inK)HO)) = i ^ (^(t + t,)(^A')t(t,)) , 

where, again, we sum the correlator over all time slices tg 
and average it. Exploiting the relation Ci2{t) = C2i(t), 
we can obtain the off-diagonal matrix element Ci2. In the 
following analysis we replace matrix element C12 with the 
matrix element C21. 

We construct the K-corrclator 

C22it)^{KHt + t,)Kit,)), (37) 

where, also, we sum the correlator over all time slices tg 
and average it. 

In this work, we also measure two-point correlator for 
the pion and kaon, namely, 

G.(t;p) = (0|7rt(O,t)7r(O,t5)|0), 

GKit;p) = (0|i^^(p,t)i^(p,ts)|0), (38) 

where the G^(t;0) is correlation function for the pion 
with zero momentum, and the G'i<-(t;p) is correlation 
function for the kaon with the momentum p. 
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IV. SIMULATION RESULTS 
A. Time correlation function 

In Fig.[3]thc individual ratios are displayed as the func- 
tions of t. The values of the direct amplitude is quite 
close to unity, indicating a weak interaction in this chan- 
nel. The crossed amplitude increases linearly, suggesting 
a repulsion in this channel. After an initial increase up to 
t ~ 4, the rectangular amplitude demonstrates a roughly 
linear decrease up to i '-^ 15, implyin g an attractive force. 
These features are what we wanted [33, S2l • 
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FIG. 3: (color online). Individual amplitude ratios it) for 
TiK four-point function calculated by wall source as functions 
of t. Direct diagram shifted by 0.8 (diamonds), crossed dia- 
gram (octagons) and rectangular (squares) diagrams. 



In Fig. m we show the real parts of the diagonal com- 
ponents {ttK — >■ -kK and k — > k) and the real part of 
the off-diagonal component ttA' — > k of the correlation 
function C{t) in Eq. (j21|) . We calculate two eigenvalues 
K{t,tR) {n = 1,2) for the matrix M{t,tj^) in Eq. ((55)) 
with the reference time t^j = 6. In Figure [5] we display 
our lattice results for A„(t, tfi){n = 1,2) as a function of 
time t. By fitting A„(f, tfl)(n = 1, 2), we then extract the 
energies which are used to calculate the scattering phase. 

We must consider the contaminations from the excited 
states and the warp-around contributions [i^l ■ By denot- 
ing a fitting range [imin,imax] and varying the values of 
the trnin and tmax, wc Can obtain the energies reliably [44 1. 
In our concrete fitting, we set tmin = tR + I and increase 
the reference time slice tn to suppress the excited state 
contaminations Moreover, we select <max to be away 
from t = T/2 to avert the warp-around effects The 
fitting parameters: tn, imin and imax are tabulated in 
Table m The fitted values for £„ (?^ = 1, 2) together with 
the fit quality x^/dof are also given in Table HI 

The energy of the free pion and kaon Ei is calculated 
from the mass TOtt and energy Ek obtained by a single 
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FIG. 4: (color online). Real part of the diagonal compo- 
nents {ttK — ^ nK and k — >■ k) and the real part of the off- 
diagonal component ttK k of the time correlation function 
C{t). Occasional points with negative central values for the 
off-diagonal component ttK — > k are not plotted. 
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FIG. 5: (color online). The eigenvalues Ai(f, t/i) and A2(t, tfl^) 
Occasional points with negative central values for the corre- 
lator \2{t,tR) are not displayed. 



TABLE I: The fitted values of the energy eigenvalues En (n — 
1,2) in lattice units. Here we also show the reference time in, 
the fitting range, tmin and tmax, the fit quality x^/dof. 













xVdof 


n 


tR 


^min 


^max 


aE„ 


1 


5 


6 


14 


0.7834(18) 


12.2/8 


2 


5 


6 


14 


0.9144(37) 


6.44/5 



exponential fit to GT,{t]0) and Gxitip) in Eq. ([55)) . as 
El = niTr + Ek , which is also listed in Table [III In Ta- 
ble [nij we show the mass m and the energy E of the pion 
and K meson with momentum p = (27r/L)e3, calculated 
from the corresponding time correlation functions [20| . 
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TABLE II: The energy of tvK system En and the s-wave 
scattering phase shift So- Ei is the energy of the free pion and 
kaon. E„ is obtained from the fitting. The invariant mass y^, 
the scattering momentum k and the phase 5q calculated with 
relations (|39II40|I in the continuum are called Cont, and those 
obtained with the relations (|41II42|I on the lattice are called 
Lat. The momentum ko is calculated by fco = 1/4 x {y/s + 
(m^ — ?7i|-)/y^)^ — m^. All values with the mass dimension 
are in lattice units. 





n = 1 




n = 2 




En 


0.8039(15) 








En 


0.7834(18) 




0.9144(37) 






Cont 


Lat 


Cont 


Lat 




0.6779(21) 


0.6888(21) 


0.8258(41) 


0.8377(41) 




0.0112(7) 


0.0150(7) 


0.0652(16) 


0.0716(18) 


kl 




0.0148(7) 




0.0700(16) 


tan (5 


0.871(34) 


0.591(22) 


-1.202(43) 


-1.495(55) 


sin^ 5 


0.431(19) 


0.259(14) 


0.591(17) 


0.691(16) 



TABLE III: Mass m of the tt, K, and k meson, 
E of the K and k meson with momentum p = 
extracted from the correlation function. 



am 
aE 



0.2459(3) 



JK 

0.3962(2) 
0.5580(14) 



and energy 
- (2^/L)e3, 



0.758(33) 
0.892(22) 



B. Lattice discretization effects 

We should consider the important discretization error 
in the t:K system's Rummukainen- Gottlieb formula (|19p. 
It stems from the Lorentz transformation from the mov- 
ing frame to the center-of-mass frame using the Lorentz 
symmetry in the continuum limit 



^MF 



p 

mz - mi. 



(39) 
(40) 



for the invariant mass ^/s, the energy in moving frame 
Emf and the momentum k. However, on the lattice, the 
discretization effects explicitly break Lorentz symmetry 
and Eqs. (p9ll40|) is only valid up to the discretization 
errors. Therefore, the definitions of ^/s and k contain 
the similar discretization errors. 

These systematic errors have been studied by Rum- 
mukainen and Gottlieb [10] , and they suggest the lattice 
modified relations. Following their suggestion, we rec- 
ommended the invariant mass ^/s and the momentum k 
from the energy in moving frame p for ■nK system as 



cosh(Vs) = cosh(£'M_F) - 2sin^(p/2) , 



(41) 



2sin^(fc/2) = cosh ( 



2^/i 



-cosh(7Ti^) ,(42) 



and then scattering phase shift was obtained by plugging 
the momentum k into the formula in Eq. (|19p. 

To understand these discretization effects, we calcu- 
late the phase shifts from the energy momentum rela- 
tions both in the continuum ([39ll40| and on the lattice 
(HUS^l)- We call the difference coming from two choices 
as the discretization error. These results are listed in Ta- 
ble |lT] along with the invariant mass yfs and the scattering 
momentum k. 



C. Extraction of resonance parameters 

From Table HTJ we note that the considerable differ- 
ences due to the choices of the energy-momentum rela- 
tions are observed in y/s and k. Moreover, the difference 
for the scattering phase 5q is significantly larger than the 
corresponding statistical errors. These are also shown in 
Fig. [6l where the phase shift sin^ 5o is also drawn. In 
Table [ll] we watch that the sign of the scattering phase 5 
at -y/s < TOk (am„ = 0.758(33)) is positive, and that at 
•\/s > m„ is negative. These features confirm a resonance 
at a mass around k mass m„. 

To avert the direct measurement of the decay 
width [U, we adopt an alternative approach. As we 
discussed in Section III Al we parameterize the resonant 
characteristic of the s-wave phase shift 5q by means of 
the effective k — >■ -kK coupling constant gKTrK, namely. 



tan<5n = 



Vi(M| 



(43) 



where Mji is the resonance mass. The equation p3|) al- 
lows us to solve for two unknown parameters: g^iTK, Mr, 
if we assume that the coupling constant g^TTK varies quite 
slowly as the quark mass changes [l^. From Table HTl we 
can notice that the difference between k and kg is not 
remarkable. In practice, we employ the scattering mo- 
mentum fco when applying Eq. (|43p . 

The lattice results of the coupling constant g^TrK and 
the resonance mass Mn solved by Eq. (|43|) are 



4.85(86) GeV, 
Mr = 0.742(21), 
Mn/m, = 0.980(51), 



(44) 

where we use the expressions p91l40p in the continuum, 
and the n meson mass is quoted from our previous 
study [lOl. If employing the expressions (|411l42p on the 
lattice, we have 



9k.itK 
Mr 

Mr/tti^ 



4.54(76) GeV, 
0.779(27) , 
1.028(58) . 



(45) 



This value of coupling constant ^kitK is in agreement with 
g^^i^ = 4.94(7) GeV obtained in Ref. [H. 

In Figure [6] we display the curves for sin^ i5o obtained 
by Eq. (j43p with the coupling constant and the res- 
onance mass Mr given in Eq. (|44p and Eq. respec- 
tively. The position at sin^ Sq — 1 which stand for Mr 
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FIG. 6: (color online). The scattering phase shift sin'^Jo, 
positions of k mass w-k and resonance mass Mr. Cont refer 
to the results achieved with the relations in the continuum 
P9ll40p and Lat to those with the relations on the lattice 
(|41ll42p . Two lines are obtained by Eq. H43p with parameters 
Qk-kK and Mr given in Eq. (|44|l and Eq. (|45|) . respectively. 
The abscissa is in lattice units. 



is also marked in Figure |6] for two cases (black cross and 
red plus for the continuum and lattice case, respectively). 
For visualized comparison, we also draw the kappa mass 
TOk (fancy cyan plus). We can note that Mfi is in rea- 
sonable agreement with the kappa mass m^. 

Assu ming that the dependence of g^TiK on quark mass 
is small [l9| , we can roughly estimate the k meson decay 
width at the physical quark mass as 



pphy 



2 

9kttK 



^phy 



Btt 



(46) 



where ml^^ = 849 MeV is the physical K meson mass, 
which we take from the most recent BESII experiment Q , 
and momentum k^^^ is calculated by 



(;jPhy)2 




where m^^^ is physical pion mass {m^^^ = 140 MeV) ij , 
and is physical kaon mass {m^^^ = 494 MeV) l|. 

This produces 



FP'^y = (335 ± 118) MeV 
where we use the data given in Eq. (|44|) . and 
rphy ^ (293 ± 101) MeV 



(47) 



(48) 



where we employ the data given in Eq. (|45p . These lattice 
estimate is fairly near to the corresponding most recent 



BESII experimental data [Bj for the n — >■ ttK decay width, 
r„ = 512 ± 80 MeV. We can observe that the difference 
stemming from our two chosen energy-momentum rela- 
tions is comparable with the statistical error. This is 
a quite inspiring result, considering that we assume the 
coupling constant is independent of the quark mass, con- 
duct a long chiral extrapolation, etc. 



CONCLUSIONS AND OUTLOOKS 



In this work, we have performed out a direct lattice 
QCD calculation of the s-wave ttK scattering phase for 
the 1=1/2 channel near the k- meson resonance region 
in the moving frame, for the MILC "medium" coarse 
{a = 0.15 fm) lattice ensemble in the presence of 2 + 1 
flavors of the Asqtad improved staggered dynamical sea 
quarks. We employed the same technique in Ref. |34ll to 
calculate all three diagrams categorized in Ref. [33, and 
obtained the pretty good signals. We have demonstrated 
that the calculation of the s-wave scattering phase shift 
for the 1 = 1/2 ttK system and then estimation of the 
decay width of k meson are practicable. The scatter- 
ing phase data shows the presence of a resonance at a 
mass around the k meson mass obtained in [2^. This 
resonance can be reasonably identified with the k meson. 
Moreover, our extracted the k meson decay is fairly com- 
pared with the K meson decay width of the most recent 
BESII experimental measurement. 

However, we realize some important issues which 
should be cleared in the future works. One is to reduce 
the discretization errors, which we show in previous sec- 
tion are comparable with statistical errors. An simple 
way to solve this problem is to use a lattice configuration 
closer to the continuum limit. Another important topics 
is to suppress the contaminations of the s-wave scattering 
phase from the p-wave scattering phase or higher, which 
we are preliminarily discussed for ttK system in Ref. [s^ . 

We adopted the effective range formula, which al- 
lows us to use the effective k — > ttK coupling con- 
stant g^TiK to extrapolate from our lattice simulation 
point (m^ -|- mK)/'mK = 0.8 to the physical point 
(to^ +mK)/mi^ = 0.73, supposing that (JkitK is indepen- 
dent of the quark mass. This is just a crude estimation, a 
more direct evaluation of the decay width is highly desir- 
able, the decay width can be estimated directly from the 
energy dependence of the phase shift data by fitting the 
BWRF if we have the simulation data which have sev- 
eral energy near the resonance mass as it was done for 
the calculation of the p meson decay width in Ref. [3] . 

Although a precise determination of the k resonance 
parameters on the lattice is absolutely a big challenge, 
our preliminary work reported here can serve as stepping 
out the solid first step in an attempt to study the strong k 
decay directly from lattice QCD in a conceptual manner. 



9 



Acknowledgments 



The definition of zcta function Zqq{s; q^) in Eq. ([8|) is 



We thank Carlcton DcTar for kindly providing us the 
MILC gauge configurations used for this work and the 
fitting software to analyze the lattice simulation data. 
We should thank Naruhito Ishizuka for the constructive 
help about group symmetry. We are indebted to MILC 
Collaboration for using the Asqtad lattice ensemble and 
MILC codes. We are grateful to Hon Qing for his sup- 
port. The computations for this work were carried out at 
AMAX, CENTOS and HP workstations in the Institute 
of Nuclear Science and Technology, Sichuan University. 



Appendix A: The calculation method of zeta 
function 



(Al) 



where the summation for r is carried out over the set 
Pd = {r|r = 7-' (n+|d) n e } , (A2) 



where a = I 



p2 



— . The operation 7 ^ is denoted 



in Eq. ((TT|) . Without loss of generosity, we consider that 
the value q^ can be a positive or negative. 

First wc consider the case of q^ > 0, and separate the 
summation in Zqq{s; q^) into two parts as 



Here we briefly discuss one useful method for the nu- 
merical evaluation of the zeta function Zoq{s; q^) defined 
in Eq. ([5]) in the center-of-mass system or Eq. ([T^ in 
the moving frame for any value of q^ (i.e., negative or 
positive). Here we follow the original derivations and 
notations in Ref. 1451. 



E 

rePd 



(r^ — q^Y 



E 



(r2 - g2).s 



E 

r^>q 



, (r^ — q^Y 



(A3) 

where the summation over r is carried out with r S P"^ in 
Eq. (jA2p . The second term can be written in an integral 
form, 



2. (,2_,2). - r(.) 2^^[j, +J, 



m Jo 



E 



3 (r2 - q'^Y ^ (^^ - 1^) 



-(r -q ) 
„2 _ „2'\s ■ 



(A4) 



The second term neatly cancels out the first term in Eq. (|A3[) . Next wc rewrite the first term in Eq. (|A4[) by the 
Poisson's rcsummation formula as 



1 



r(5) Jo 



(A5) 



The divergence at s = 1 comes from the n = part of the integrand. Therefore we divide the integrand into a 
divergent part (n = 0) and a finite part (n ^ 0). The divergent part can be evaluated for Res > 3/2 as 



3/2 ^21 



3/2 n 



(A6) 



The right hand side of this equation has a finite value at s = 1. 

After connecting all terms we obtain the representation of the zeta function at s = 1, 



VI^-Zi,il;q^)^J2- 



j,2 _ q2 



At e*«" 



7r\ 2 
1 



/=0 



/2 l\ 



(A7) 



where X^Lgz^ stands for a summation without n = 0. 

For the case of q^ < 0, it is not necessary to divide the summation in Zoq{s; q^), and it can be also expressed in an 
integral form. 



E 



E 



iri - q^Y ^ (r2 - q^Y r(s) 



(A8) 
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Following the same procedures, we arrive at the same expression in Eq. (|A7[) . Hence, Eq. (jA7[) can be applied for 
both cases. 

Substituting d = (1,0,0) into Eq. ()A7|) we obtain the representation of the zeta function appeared in Eq. (f20|) 



where the imaginary parts are neatly canceled out. 

Substituting d = and 7 = 1 into into Eq. (|A7|) . we obtain the representation of the zeta function in the CM 
system appeared in Eq. ([5]) 

^M-E'-^- [ (7)' + £ ^ (Aio) 

I 



I also note that, for negative g^, an asymptotic expres- of the zeta functions with that of above described repre- 
sion of the zeta function 2^oo(s;?^) has been derived in sentation, and found reasonable agreement. 
Rcf. [31[. We numerically compared this representation 
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